Origin of the efficient light emission from inversion domain boundaries in GaN 
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Intentionally-produced inversion domain boundaries in GaN have been reported to be highly 
efficient shallow recombination centers. Here I report a rationale for this phenomenon based on ab 
initio density-functional calculations. I also propose a model, based on the existence of polarization 
in GaN, of the observation that a domain boundary acts as a rectifying junction under voltage 
applied between the two opposite-polarity surfaces. 

PACS numbers: 61.72.Bb, 61.72.nn, 77.22.Ej 



Nitride semiconductors, and GaN in particular, are 
major competitors in the optoelectronics market at the 
high frequency end of the visible spectrum. Understand- 
ably, considerable attention has been devoted to the na- 
ture of recombination centers in In-alloyed GaN, both 
non-radiative (mainly threading dislocations) and radia- 
tive (allegedly. In-rich regions). Not nearly as much work 
exists on recombination centers in plain GaN. In recent 
times, extended defects originating from the formation of 
polarity domain boundaries, and known as inversion do- 
main boundaries (henceforth IDB for brevity) , have come 
into focus in this context. They have been studied theo- 
retically directly imaged by high-resolution TEM ||] 
and piezoresponse force microscopy , and finally, inten- 
tionally produced in MBE epitaxial growth on patterned 
substrates . By means of integrated PL topography Q 
it was observed that, unexpectedly, these IDBs exhibit a 
tenfold more intense light emission compared to plain 
Ga-face surface areas. The transition energy is about 30 
meV below the normal emission, and unlike the latter, 
it is polarized in the plane of the defect. The IDBs thus 
qualify as highly efficient recombination centers. The 
first purpose of this paper is to report a rationale for this 
phenomenon, based on ab initio density-functional calcu- 
lations. A further puzzling observation is that an IDB 
between metal contacts placed on the Ga-face and N-face 
surface regions, acts as a rectifying junction (direct bias 
being Ga-face to N-face) . Here, I also propose a model of 
the latter behavior, based on the existence of permanent 
bulk polarization in GaN ||]. 

Being a wurtzite crystal, GaN possesses a singular po- 
lar axis labeled (0001) with the convention that the pos- 
itive direction points from the cation to the anion, and 
an attendant permanent polarization vector oriented as 
(0001). There exist two geometrically distinct crystal ter- 
minations, the (0001) or Ga-face and (0001) or N-face, 
whereby a Ga (respectively, a N) atom would have a sin- 
gle "dangling bond" orthogonal to the surface. Two re- 
gions of opposite polarity, conveniently labeled Ga-face 
and N-face, meet at an IDB. There are at least two 
plausible models thereof, which were discussed earlier by 
Northrup et al. [0. In the first, cations and anions are 



exchanged in half of the crystal, giving rise to a geomet- 
rically undistorted interface, which however contains two 
(one N-N and one Ga-Ga) wrong bonds per ceU §. The 
second model, also known as IDB* and found in Ref.|l| 
to be decidedly favored energetically over the first, is ob- 
tained shifting one of the inverted half-cristals by half the 
cell length (c/2) along the singular axis. This enables all 
atoms to be correctly bonded, at the only cost of a mod- 
est distortion of the local bond angles, the largest single 
angle change being ^^109° to 90°. In this work, I will be 
considering the latter model. 

I calculate total energies, forces, stress, and elec- 
tronic potential and density from first-principles, within 
density functional theory in the Ceperley-Alder j|] 
local-density approximation to exchange and correlation, 
using the accurate all-electron frozen-core Projector- 
Augmented Waves (PAW) method as implemented 
1 11 1 in the Vienna Ab-initio Simulation Package (VASP) 
|12|. The plane wave cutoff is set at 290 eV. The par- 
tial wave expansions for N and Ga include angular mo- 
menta up to 1=2 and /=3, respectively, with all on- 
site projectors included. Ga Sd electrons are treated 
as valence. Accurate reciprocal-lattice grids are used, 
with no wrap-around aliasing allowed in Fourier trans- 
forms. Monkhorst-Pack k-space integration grids are 
used; for the bulk, converged results are obtained with 
the (4x4x4) mesh, tested up to the (12x12x12). There- 
suiting crystal structure of GaN is aGaN=3.160 A, c/acaN 
=1.626, u=0.3769, which matches closely (maximum de- 
viation: 0.9 % for flGaN) the accepted experimental val- 
ues acxp=3.189 A, c/acxp=1.626, u=0.377 Q - one of 
the best descriptions of bulk structure ever achieved for 
GaN. 

The IDB* planar defect can be simulated in periodic 
boundary conditions via repeated supercells. I use a crys- 
talline GaN orthorombic supercell comprised of 24 for- 
mula units (48 atoms), with linear dimensions (1, 1.626, 
6 y/i) in units of acaN- After testing the (nxlxn) and 
{nx2xn) k-point meshes with n=2,4, and 8, I used the 
(4x2x4) for all the results reported here. The symme- 
try of the system is found to be almost exactly T)2h^ al- 
though this symmetry is not explicitly enforced. When it 
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is, the k-space grid includes 4 special points in the super- 
cell Brillouin reduced zone; otherwise the mesh comprises 
8 points in the reduced zone. 

Before relaxation, the formation energy of the IDE de- 
fect is 40 meV/A^. To allow for possible defect-induced 
strains, the cell length in the y direction (orthogonal to 
the defect plane) was adjusted, with the atomic geome- 
tries being constantly optimized to obtain residual force 
components below 0.01 eV/A. The cell length does not 
change within numerical accuracy; some local relaxation 
occurs in the planes immediately adjacent the defect, all 
other atomic positions being unperturbed. Relaxation 
reduces the formation energy to 20 meV/A^, a rather 
small value compared to the 114 meV/A^ of two relaxed 
(lOTO) surfaces of GaN |l|l, and agreeing well with the 
25 meV/A2 of Ref.§. 

I calculate the transition energies, using the so called 
ASCF principle [Q, as the combination of (variational) 
total energies of the system with varying number of elec- 
trons 

i^gap = EtotiN -hi) - 2Et,t{N) + Et,t{N - 1). (1) 

This expression is exact in the N ^ oo limit, but it turns 
out to be pretty accurate also at finite TV to estimate 
the gap of crystals jl^. The added charge is compen- 
sated by a uniform background. Since part of the excess 
charge localizes in the defect plane, I apply a correction 
CT^ c?ccii/2er the interface charge, dccii the cell length, 
Er the static dielectric permittivity) to account for the 
spurious interaction of such sheet charges, screened by 
the static dielectric response of GaN. This correction is 
in the order of 10 /xeV, i.e. fully negligible on the scale 
of other computational approximations (e.g. frozen core, 
LDA, etc.) and convergence errors. 
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FIG. 1: Local self-consistent electronic potential at the IDB* 

The transition energy obtained for the bulk undefected 
cell is 3.49 eV, in good agreement with experiment. For 



the defected cell, the transition energy is also 3.49 eV, 
thus predicting a red shift of zero, compared to the ex- 
perimental value of 30 meV. In the calculations the excess 
electron or hole are added separately to the system, and 
they do not interact in any way, so that excitonic effect 
are not described. 

However, in view of the low experimental temperature, 
it is fairly natural to suppose that the observed shift may 
be of excitonic origin. Although electron-hole interaction 
is not included, it is possible to pinpoint the origin of the 
effect looking at the precursors of the exciton, i.e. the 
electron and hole densities and attendant self-consistent 
potential. Specifically, let us examine in Fig. |] the self- 
consistent electronic potential of the defected cell in the 
neutral charge state (no extra charge), filtered through 
a macroscopic average [p^ eliminating the periodic mi- 
croscopic oscillations. As can be seen, the potential in 
the vicinity of the defect is such as to attract electrons 
with its weak but far-reaching side wings on the outer 
side of the defect core, and repelling them off the core it- 
self; contrarywise, holes are attracted to the defect core. 
Both electrons and holes are thus expected to localize in 
the vicinity of the defect (which should enhance excitonic 
interactions). 




FIG. 2: Density of excess hole (left) and electron near the 
IDB* defect. 

This is confirmed by inspection of the excess electron 
and hole charge densities, obtained as density differences 
between charged and neutral cells. As can be seen |l^ in 
Fig. ^, electrons accumulate just outside the defect and 
onto the central repulsive wall, whereas (although this is 
not especially clear from the picture) the hole density is 
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enhanced in the defect core. 

Inspection also helps rationalize the observed emission 
polarization in the defect plane. The electron density, 
mostly of cation s character, is "planarized" within the 
two-dimensional lateral wells of the defect potential. The 
hole density away from the defect consists of a lobe along 
the y direction: this is because the valence top state in 
GaN is a Fg doublet of PxPy nature (the x lobe is not 
visible, as the x direction is orthogonal to the plane of the 
picture) ; near the defect the py lobe pivots around the x 
axis admixing some character of the pz singlet Fg which 
lies some 50 meV lower than Fg in unperturbed GaN. 
Overall this results in an enhanced directional overlap of 
the electron and hole functions. 

A rationale for this admixture can be extracted from 
the following model (refer to Fig. ^ . Consider zincblende 
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FIG. 3: Schematic of valence states admixture between 
wurtzite-like and NaCl-like environments. 

GaN, with its three-fold degenerate valence band top. 
Stretch a bond along one of the (lll)-equivalent direc- 
tion; the degeneracy is removed, and as in wurtzite the Fg 
singlet is lower than the Fg doublet. Now displace, say, 
the anion along the chosen (111) direction, and monitor 
the relative position of the states. What is observed is 
the situation schematized in Fig. |[ as the distortion pro- 
gresses the Pz singlet nears the doublet, becomes higher 
in energy (i.e. the band top), then merges back with 
the doublet into the threefold-degenerate valence top of 
the NaCl structure. The IDB* embodies indeed a lo- 
cal NaCl-like geometry, with one of the a-plane bonds at 
right angles with the polar-axis bond. This distortion is 
the cause of the admixture of the two top valence states. 
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Another remarkable effect, observed by Stutzmann et 
al. 1^, is that an IDB* acts as a rectifying junction 
when electrically biased by two contacts placed respec- 
tively on the adjacent Ga- and N-face regions, direct- 
bias being realized by the positive pole on the Ga-face 
region. This behavior can be qualitatively explained by 
the existence of an intrinsic macroscopic polarization in 
wurtzite as schematized in Fig. ^. The two adjacent 
opposite-polarity regions possess opposite intrinsic polar- 
ization vectors. This results in fixed polarization charges 
at the upper and lower ends of the sample, and chang- 
ing sign across the IDB junction. On the upper side, 
mobile charges or impurities from the vacuum neutralize 
the polarization charges on each of the regions. At the 
border between the latter, the compensating charges will 
mutually neutralize and deplete, creating a junction-like 
potential which opposes the current flow from Ga-face 
to N-face. Direct polarization with the positive pole on 
the Ga-face side will then cause a current flow from Ga- 
face to N-face. It is difficult to say whether this current 
will flow subsurface or otherwise. Certainly, however, 
the lower side of the polarized layer, a few /im below the 
surface, will not be involved. Thereby, charges from, or 
defects in the nucleation layer are assumed to neutralize 
permanently the polarization charge. 

In summary I have given an interpretation of the ob- 
served enhanced light emission from inversion domain 
boundaries in GaN, and I have suggested a qualitative 
explanation of the rectification effect observed across the 
boundaries themselves. 

I dedicate this work to the memory of my father Ugo 
Fiorentini. Support from the Italian Ministry of Re- 
search under the PRIN funding programme (years 2000 
and 2001) is acknowledged. 



FIG. 4: Schematic of the polarization-driven rectifying mech- 
anism. 
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